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We study the dynamics of an impurity of finite mass injected into a one-dimensional quantum 
liquid at zero temperature, either at finite velocity or at zero velocity with a force driving the 
impurity. We obtain accurate results using numerical simulations based on matrix product states, 
and find that in both cases, the impurity undergoes oscillations, however the physical mechanism 
is different: the driven impurity undergoes Bloch oscillations by following the ground state branch 
while continuously emitting phonons, whereas the undriven impurity undergoes coherent quantum 
oscillations at an emergent energy scale, called quantum flutter in previous work, whose amplitude 
grows with increasing initial velocity. We find these results to be independent of whether the system 
is integrable or not, and robust to changes in the microscopies of the model, suggesting that they 
are universal. 



I. INTRODUCTION 

Experiments in ultracold atomic systems have recently 
realized different incarnations of a quantum impurity 
problem in which, starting from a ID gas of particles 
with the same internal state called the background gas, 
a single particle called the impurity is promoted to a 
different internal state and its properties over time are 
studied [1—5]. The background itself exhibits behavior 
special to one-dimensional quantum systems [6, 7]. In- 
vestigations of mobile impurities in such systems con- 
tributed to the understanding of various phenomena in 
ID many-body physics: excitation spectrum and effective 
mass [8-11], orthogonality catastrophe [12, 13], logarith- 
mic diffusion of Green's functions [14, 15], persistence of 
threshold singularity in spectral functions [16, 17] and its 
momentum-dependent power-law scaling [14, 18-22], and 
response to an external driving force [23-27] and external 
confinement [28-30]. 

Recently it was shown that if an impurity, initialized 
in a plane wave with a finite momentum Q, is added to 
a gas of free fermions (or bosons interacting with each 
other through infinitely strong (^-function potential) the 
impurity only sheds a part of its momentum to the back- 
ground gas and forms a correlated state that no longer 
decays [31]. Furthermore, if Q is of the order of or larger 
than the Fermi momentum kp the momentum of the im- 
purity undergoes coherent oscillations which were called 
quantum flutter and claimed to be due to the impurity 
forming a coherent superposition of families of quasi- 
equilibrium states called polaron-like and exciton-like. 
These phenomena were demonstrated by examining the 
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FIG. 1. Schematic picture of the setup: (Color on- 
line) An impurity atom, blue sphere, moves with momentum 
{P±(t)) through an interacting bosonic background gas ini- 
tially in its ground state, red spheres. The impurity is in- 
jected either in a plane wave with momentum (Pj.(0)) = Q 
or at zero momentum with a constant force of strength e a 
driving it. The impurity interacts with the background par- 
ticles through a contact interaction of strength 7, and the 
background particles interact with each other with a contact 
interaction of strength -jbg ■ 



full quantum-mechanical evolution obtained from the ex- 
act Bethe Ansatz solution which exploits the integrability 
of the system. Integrability implies the existence of an 
extensive number of mutually commuting integrals of mo- 
tion [32, 33]. Thus integrability may largely constraint 
the dynamics of a system, which poses the question how 
universal the results obtained for exactly solvable mod- 
els are. As a general rule, the low-energy dynamics of ID 
gapless quantum systems is universal and does not differ- 
entiate integrable and non-integrable systems [6, 7], while 
the equilibration after a strong perturbation reveals a 
vast amount of integrability-specific phenomena [34-36] . 
Stemming from a far-out-of equilibrium initial configura- 
tion, quantum flutter is a natural candidate to be tested 
against various integrability breaking conditions, as we 
do in this work. 

Since the undriven impurity undergoes dissipationless 
motion, one question which naturally arises is what the 
effect is of a driving force acting only on the impurity. 
A semiclassical analysis [27] predicted that the system 
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FIG. 2. Quantum nutter: (Color online) An impurity 
(down spin) injected in a plane wave with momentum Q into 
the background gas (up spins) initially slows down, however, 
it never comes to a full stop. Further, the impurity mo- 
mentum {P_i(t)) undergoes oscillations in time with an am- 
plitude that grows with increasing initial momentum Q — a 
phenomenon known as quantum flutter, whose full solution for 
the integrable case of an infinitely repulsive background gas 
has been explored in Ref. [31]. In this work, we demonstrate 
that quantum flutter is not specific to integrable systems and 
robust to changes in the microscopic model. As an example 
this figure shows quantum flutter in the impurity momentum 
for a background gas with finite interactions 7b 9 thus break- 
ing integrability. Units: the background gas has density p-f 
which defines the Fermi momentum k F = Trpf Momentum 
is given in units of fc^, and times in units of the Fermi time 
tF — 2m^/k F where is the mass of the background parti- 
cles. Here the masses are considered to be equal: m-f = m\.. 

would reach a finite drift velocity and undergo Bloch 
oscillations around this velocity. Intuitively, the phys- 
ical picture is that the impurity is constantly emitting 
phonons in the background gas and following the ground 
state branch, which is generically periodic in momentum 
space for one-dimensional systems, thus leading to peri- 
odic effects in the dynamics. The semiclassical analysis 
assumes that the impurity has a well defined position as a 
function of time, which implies that the wavefunction of 
the system can be written as a product of wavefunctions 
for the impurity and for the background gas. This ansatz 
should be accurate for heavy impurities, however, for im- 
purities of mass around that of the background particles 
assessing the validity of that assumption requires an ac- 
curate simulation of the time evolution of the system. 

In the present paper we present essentially exact re- 
sults for the dynamics of an impurity injected into a ID 
quantum gas by employing the time-evolving block deci- 
mation (TEBD) algorithm [37, 38]. First, we show that 
quantum flutter is not restricted to integrable systems 
and is robust to changes in the microscopies of the model 
(see Fig. 2 for representative results). Second, we inves- 
tigate the effects of a finite force acting on the impurity 
and find oscillations in the impurity momentum which 
decay in time and oscillate at a frequency higher than 




FIG. 3. Bloch oscillations: (Color online) (a) Schematic 
dispersion relation of a repulsive impurity in a background 
gas, lower solid curve, above which a continuum of manybody 
excitations exists, that are composed of the moving impurity 
dressed by a certain number of phonons. For comparison the 
schematic dispersion of the background gas is shown as well, 
upper dashed curve, which is linear at low momentum. A 
force, which acts adiabatically on the impurity, would drive 
it through the ground state branch leading to Bloch oscilla- 
tions with period 2k F without energy transfer to the back- 
ground gas. In contrary, for a finite force the acceleration of 
the impurity is accompanied with the radiation of phonons, 
which effectively heats up the system and gives rise to a fi- 
nite drift velocity around which the impurity performs Bloch 
oscillations [27]. (b) shows exemplary curves for the impu- 
rity momentum as a function of time for two different values 
of the force e a confirming the picture of (a). Furthermore, 
the amplitude of Bloch oscillations around the drift velocity 
is damped due to nonlinear effects. A stronger force e a gives 
rise to an increasing drift velocity va as well as an increasing 
oscillation frequency w&, as shown in the inset. For details see 
main text. 



predicted on the basis of a semiclassical analysis [27], 
see Fig. 3. We compare our results with the experimen- 
tal results from Ref. [1] which explores such a situation 
and find that they are quantitatively consistent. We also 
explore the evolution of the background density around 
the impurity and find signatures of the driven impurity 
problem which could also be studied in experiments. 

This paper is organized as follows: In Sec. II we present 
the model at study and discuss details on the numerical 
method. Results in the absence of a driving force are 
discussed in Sec. Ill, where we show that experimentally 
relevant integrability breaking terms do not affect the 
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qualitative features of quantum flutter. In Sec. IV we 
discuss the emergence of Bloch oscillations for an im- 
purity which is driven by a constant force. Finally, we 
conclude our findings in Sec. V. 



II. THE MODEL 

In the model we consider a quantum particle (im- 
purity) embedded into a gas of quantum particles of 
other type (background gas), see schematic illustration in 
Fig. 1. The motion of the impurity and of the background 
gas is confined to one spatial dimension. We perform the 
numerical calculations using a lattice model with a finite 
number of sites N s . The system is always studied at low 
filling, so that the obtained results are representative for 
the continuum. The Hamiltonian of the model reads 



We denote the ground state of the background gas, 
Eq. (2), as |BG). We then initialize the impurity in a 
plane wave with momentum Q: 



m = 0))=J Qi \BG), 



(6) 



where the operator Cq^ = 
impurity with momentum Q. 
tions the state \tp(t)) = e~ lHt 
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where 
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describes the background gas, 
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describes the impurity, and 
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Ylj=i e iC ^ Xs Cj, creates an 
In our numerical simula- 
\ip(t — 0)) is represented 
as a matrix product state (MPS). We start with deter- 
mining the initial ground state |BG) by density matrix 
renormalization group [39, 40], and then add the impu- 
rity as given by Eq. (6). Subsequently, we evolve the 
system in time using the time-evolving block decimation 
(TEBD) algorithm [37, 38]. We push TEBD to perform 
high-accuracy simulations: Specifically, the presented re- 
sults are obtained for systems with N s = 400 lattice sites, 
N = 40 particles, and MPS bond dimension M = 400 or 
M = 500, which determines the numerical accuracy of 
the calculations. In Fig. 4, we consider even a double as 
large MPS bond dimension M = 800, which increases the 
numerical effort by a factor of eight, and in Fig. 5 we in- 
creased the simulated systems to N s = 600, N = 60, and 
M = 600. We verified that all the results are representa- 
tive for the continuum and when appropriately rescaled, 
see below, do not depend on N, N s and on the MPS 
bond dimension M. 

We are aiming our numerics at the analysis of the con- 
tinuum limit of the model defined by Eq. (l)-(5). The 
limit is taken as usual (see, for example, Chapter 2 of 
Ref. [41]). Namely, Cj a — > y/a^ a (x), where x = aj. The 
density operator in the continuum p G {x) = ^\{x)^ a {x). 
The summation over lattice sites is replaced with the 
integration, J2f=i a f{ a j) §q dx f (x), where / is an 



describes the interaction of the impurity and background arbitrary function, and L — N s a. Thus, 



gas. Here c\ al c ia , and p ix 



c\ a c ia are the creation, an- 
nihilation, and the particle number operators on the zth 
lattice site, respectively. The subscript a =t,l; we call 
the impurity a spin-down particle and the particles in the 
background gas spin-up particles. There is only one im- 
purity in the system, and the number of the background 
gas particles, N, is arbitrary, 



where m a = l/(2J a a 2 ) and 
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The parameter Uff (U^) characterizes the on-site re- 
pulsion energy between two background particles (back- 
ground particle and impurity); J„ specifies the hopping 
energy. The last term on the right hand side of Eq. (3) de- 
scribes a constant force F applied to the impurity, Xi — ia 
being the position of ith site and a the lattice spacing. 
The background particles satisfy bosonic statistics: the 
operators cL and have canonical commutation rela- 
tions, 



The impurity-background interaction term, Eq. (4) , takes 
the form 



Htf=9 / dx p l (x)p t (x), g = aU lt . (9) 
Jo 

It is convenient to use the dimensionless parameter 7 for 
the strength of the impurity-background interaction: 



-it 
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] = 5ij. Thus, the model (1) is a partic- 
ular case of the spin- 1/2 Bose-Hubbard model (Eq. (5) 
restricts the number of spin-down particles to one). 
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Here we used that p^ = N/L and = N/N s . The 
continuum limit of the interaction term in Eq. (2) is 
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FIG. 4. Impurity momentum as a function of time: (Color online) An impurity with momentum Q is immersed in a 
one-dimensional interacting bosonic liquid. The expected impurity momentum (Pj.(£)} is plotted as a function of time. {Pi(t)) 
first decays to finite value, and then oscillates as a function of time. The simulations are carried out using TEBD on a spin-1/2 
Bose-Hubbard model that is dilute enough to obtain results that are valid for the continuum limit. The interaction between 
the impurity and background particles is fixed at 7 = 12. The simulations have converged and do not depend on number of 
lattice sites N s , number of particles N or MPS bond dimension M (the largest system we simulated had N s = 400, N = 40 
and M = 800. N/N a must remain small enough to avoid lattice effects, we found the results to be independent of N as long 
as N/N s < 0.2). (a) Time evolution of (P±(t)) for a set of different initial momenta and background interaction strengths 
jb g . After the initial decay, the plots depend weakly on the initial momentum Q, provided Q is large compared to kp. Inset: 
quantum flutter amplitude A as a function of initial momentum, where the amplitude is defined as the difference in (P±(t)) 
between the first maximum and first minimum. The amplitude grows with Q, and the second derivative of A as a function 
of Q changes sign at Q = kp. The growth in A is increasingly sharp as 7^ increases, in other words as the system becomes 
more fermionic (at jb g = 00 the system is fermionizable). Colors indicate yt g as in the main figure, (b) A closer look at the 
dependence on 7& g . Plotted is (P±(t)) for 7 = 12 and for 7^ = 4, 8, 12, 20, 00 from bottom to top. The limits jt a = 00 and 
7b g = 7 are both integrable, the former via fermionization, the latter being the Gaudin-Yang model. Inset: amplitude as a 
function of fbg- The solid line, labeled Ai, is defined as in the inset of (a), and initially decays with decreasing 755 but finally 
increases again for small enough yi,g, where one sees a pronounced revival of momentum which then reaches a maximum and 
finally decays to zero while oscillating with negligible amplitude. The first minimum and maximum are therefore special, as 
the system is still forming the correlation hole around the repulsively interacting impurity and only settles into quantum flutter 
afterwards. Thus a possibly more meaningful definition of quantum flutter amplitude is given by the dashed line, labeled A2, 
which shows the difference in (Pi(t)) between the second minimum and second maximum. This quantity strictly decreases with 
decreasing 7f, 9 . This indicates that the oscillations are strongest in the fermionic limit. See also the main panel for a graphical 
illustration of the definition of the amplitudes A\ and A-z. 



=rr / dx pi(x)p-f(x) and the dimensionless interaction 
strength between the background particles is 
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The model we consider is exactly solvable for some 
particular values of parameters (as mentioned in the In- 
troduction, we interpret the model to be integrable for 
those values of parameters): A necessary condition for 
the integrability is = Jj,. When F ^ the model (1)- 
(5) (and its continuum limit) admits exact solution for 
= = +00 only. Letting F = makes the 

model (l)-(5) exactly solvable for = 00 and arbitrary 
U±f, where it can be mapped onto a model of fermions 
interacting with a single particle of other type. That 
mapping was used, in particular, in Ref. [31]. When 
U-ff = 7^ 00 the model (l)-(5) is not integrable, how- 



ever, its continuum limit is integrable: it is known as the 
Gaudin-Yang model [42, 43]. For any other choices of the 
parameters the model is, most probably, not integrable. 



III. DYNAMICS IN THE ABSENCE OF A 
DRIVING FORCE 



We consider the case without a force F = and the 
impurity initialized in a plane wave with finite momen- 
tum Q on top of the background liquid in Fig. 4. We 
assume equal masses = m-j-, unless mentioned other- 
wise. We track the expected impurity momentum (Pj,(i)) 
as a function of time, where Pi = J2k^ c ki c ki- ^ e de- 



tine the operator Tj = ct +1 . c^, — c^c i+1 ^. If the system 
were translationally invariant, one could use the fact that 
Tj = — 2i sin(fca)cL c fe , , which in the limit a — > 
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becomes —2iaP^, to calculate the impurity momentum. 
In our simulations open boundary conditions are used, 
which arc not translationally invariant. However, near 
the center of the trap (ip(t)\Tj\ip(t)) is weakly dependent 
on the site j at times before finite size effects become 
important, which occurs once excitations start to reflect 
off the boundary and arrive at the center of the system. 
We average (ip(t)\fj\ip(t)} over the central 20% of the 
simulated system, and verified that the results did not 
change with modest change of the size of the area that 
is being averaged over. From this quantity we can di- 
rectlyuse (P ; (t)> = ^ E^fig/aM^M*)) 
A = 0.2jV s to extract the impurity momentum. In gen- 
eral we find that the momentum does not decay to zero, 
instead it decays to a finite value. Furthermore, the mo- 
mentum shows damped oscillations as a function of time 
whose amplitude grows with increasing initial momentum 
Q, see inset in Fig. 4 (a), where the amplitude is defined 
as the difference in (P^(t)) between the first minimum 
and first maximum. The second derivative of the ampli- 
tude as a function of Q is zero at Q = kp . In Fig. 4 (b) 
we study in more detail the dependence of (P\,(t)) on 
the background gas interaction strength j^g- Through- 
out this work we use the Fermi momentum Uf = "Kp-\ and 
the Fermi time tp — 2m+/kp as proper units. 

In the case of jbg = oo it was argued that quantum 
flutter results from the impurity forming a superposition 
of two families of quasi-equilibrium states [31], one with 
the total system (background gas plus impurity) in its 
lowest energy state at momentum k = kp and the other 
with the dressed impurity at zero momentum k = and 
an excitation of the background at k = kp. If the ini- 
tial momentum Q > kp, any momentum in excess of 
kp is sent out in a wave packet. So far the energies of 
these states have not been calculated except in the limits 
where the system is integrable [22, 31]. Since in passing 
we obtain results for the symmetric Gaudin-Yang case 
7 = 7b 9 in this work, we verify that this relationship 
holds in that limit as well, as we plot in Fig. 5. For the 
Gaudin-Yang model the ground state of the background 
plus impurity at finite momentum is known as a magnon, 
and finite momentum excitations of the background are 
known as plasmons [14, 22]. Quantum flutter for this 
bosonic system is due to the dynamical formation of a 
coherent superposition of a plasmon at kp and a magnon 
at kp , as we describe in more detail below. 

The red squares in Fig. 5 (a) denote the frequency of 
the oscillations in (_Pj_(t)), which we call the quantum 
flutter frequency hwf, extracted from the TEBD simula- 
tions. Specifically, we take the set of distances in time 
from one extremum in (P|_(f)) to the next, which each 
give a value St for half of the flutter period. The mean 
of this dataset (St) is related to the flutter frequency 
through U3f = n / (St) shown as red squares and the corre- 
sponding error bars are determined from the propagated 
standard deviation. 

The plasmon-magnon gap at kp is indicated by the 
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FIG. 5. Quantum nutter frequency: (Color online) (a) 
Comparison of the quantum flutter frequency of oscillations 
in (P_i(t)) with the plasmon-magnon gap for the symmetric 
model, where 7 = 7b 9 . The red squares denote the quan- 
tum flutter frequency ujf extracted from the TEBD simula- 
tions while the solid line corresponds to the plasmon-magnon 
gap at kp obtained exactly with Bethe Ansatz [22], which 
the data agrees with within the error bars. The inset shows 
the plasmon and magnon dispersions at 7 = jtg = 10. The 
magnon dispersion is the ground state energy of the system 
(background gas plus impurity) as a function of total momen- 
tum k. The plasmon dispersion gives the minimum energy of 
particle-hole excitations with total momentum k in the back- 
ground gas, which is linear for small k [44]. (b) Intuition 
behind quantum flutter: An impurity injected into the back- 
ground gas with momentum Q = kF, is in a superposition 
of a magnon (a dressed impurity) with momentum k = kp, 
left, and a plasmon (a density excitation of the background 
gas) with momentum k = k F while the dressed impurity is at 
rest k = 0, right. Any initial momentum Q > k F in excess 
of k F would be emitted in a wave packet excitation of the 
background gas. 



solid curve which agrees with the flutter frequency to 
within error bars. Plasmon E p (k) and magnon E m (k) 
dispersions are shown in the inset for 7 = 7^ = 10. The 
magnon dispersion E m (k) gives the ground state energy 
of the system (background gas plus impurity) as a func- 
tion of total momentum k with E m (k = 0) set to zero. 
The plasmon dispersion E p (k) is linear at small momen- 
tum [44] and characterizes the lowest energy to create 
particle-hole excitations in the background gas at finite 
total momentum k, once again with E p (k = 0) set to zero. 
The difference between the plasmon and magnon disper- 
sions at k = kp is what we call the plasmon-magnon 
gap. 

The interpretation is as follows: if the impurity is in- 
jected in the system at momentum Q = kp, it can either 
form a magnon (a dressed impurity) at momentum kp 
or create a dressed impurity with zero momentum and 
emit a plasmon (an excitation of the background) with 
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FIG. 6. Formation of the correlation hole and wave packet emission: (Color online) Plot of the density-density 
correlator G^(x,t) = {pi(0,i)pf(x,t)), which corresponds to the density of the background gas around the impurity. Shortly 
after injecting the impurity a pronounced correlation hole (blue valley) is formed due to the repulsive interactions. This process 
is accompanied by the emission of a wave packet (red peak) in the background gas that carries away the impurity impurity 
momentum. The time evolution of the correlations Gif(x, t) is shown for 7 = 12, Q = 1.5/cf, and (a) 7^ = 00, (b) jb g = 20, 
and (c) 7(, 9 = 5. Other system parameters are the same as in Fig. 4. The qualitative features do not depend on the background 
interaction. The main changes as 7^ is lowered are the disappearance of the Friedel oscillations in the background, and growth 
in the size of the correlation hole. 



momentum kp. In Fig. 5 (b) we show the resonant case, 
where the incoming momentum is precisely enough to 
create a deep hole. If the impurity is initialized with a 
momentum Q larger than kp the excess in momentum 
is carried away through an emitted wave packet. In the 
resonant case the total energy of the magnon-like state 
is + E m (k = 0) + E m (kp) and of the plasmon like 
state E N ^ + E m {k = 0) + E p (kp ), where E N ^ is the total 
energy of the N background gas particles. Thus the dif- 
ference between E\ and E2 is the plasmon-magnon gap 
Ep(kp) — E m (kp) which corresponds to the flutter fre- 
quency. This strongly suggests that the physics is gov- 
erned by the formation of a coherent superposition of 
these two quasi-equilibrium states. 
The correlation function 



(12) 



measures the background density with respect to the im- 
purity position as a function of time, and shows a pro- 
nounced peak which corresponds to the emission of a 
wave packet that carries away most but not all of the 
impurity momentum, see Fig. 6. Furthermore the repul- 
sively interacting impurity builds up a correlation hole 
around itself. These qualitative results do not change 
as one varies the background interaction strength. The 
main changes as one goes away from the Tonks- Girardeau 
(TG) limit -fb g = 00 are that the Friedel oscillations dis- 
appear and that the width of the correlation hole in- 
creases. Note that whereas in [31] Friedel oscillations ap- 
peared after a finite time, as a consequence of the formed 
correlation hole which breaks translational invariance, 
their presence here at all times is quite natural, because of 
open boundary conditions used in the numerical simula- 
tions. In particular, when moving away from the bound- 
ary, Friedel oscillations decay as cos(kpx)x~ 2K , where K 
is the Luttinger parameter [6, 7], giving rise to a faster 



decay with increasing K, which is in agreement with the 
results in Fig. 6. 

Changing the mass of the impurity with respect to 
the mass of the background particles, which makes the 
system nonintegrable, does not affect the qualitative re- 
sults, see Fig. 7. We find that the impurity momentum 
(Pj,(i)) exhibits quantum flutter also in that limit pro- 
vided the ratio of the masses does not get too large or too 
small. Empirically, we find that for 0.5 < m^/m^ < 2.0 
quantum flutter is still pronounced, see Fig. 7. In the 
inset of this figure we show the amplitude of the oscilla- 
tions where, as in Fig. 4 (b) the solid line corresponds to 
the amplitude of the first revival of whereas the 

dashed line indicates the difference in (P^(t)) between 
the second maximum and second minimum. From the 
latter it turns out that the integrable point — ex- 
hibits the strongest quantum flutter oscillations. This is 
in contrast to the first revival, whose amplitude is shifted 
toward lighter impurities, which can be attributed to the 
formation of the correlation hole at short times. 

An experimental realization in a cold atomic setup will 
typically have a harmonic background trapping potential, 
thus we add to the Hamiltonian: 



<J 1 — 1 



(13) 



where lot is the trapping frequency. The trap breaks 
translational symmetry, thus momentum is no longer a 
good quantum number, which begs the question how 
quantum flutter is affected. In Fig. 8 the impurity mo- 
mentum is shown in the presence of a trapping potential. 
The frequency of quantum flutter ujf increases quadrat- 
ically with the trap frequency u>t, see inset of Fig. 8. In 
this inset, we plot a gray line that corresponds to an 
equilibrium calculation of an energy which agrees exactly 
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FIG. 7. Effect of mass imbalance: (Color online) When 
the masses mi 7^ uif, the model is no longer integrable. 
This does not change the phenomenology of quantum flutter. 
In this plot we take the background to be in the TG limit 
(jb g = 00), the impurity-background coupling 7 = 12, and 
the initial impurity momentum Q = 1.2&F- Other parame- 
ters are as in Fig. 4. Inset: The lines are defined as in the 
inset of Fig. 4 (b). The amplitude A? of the quantum flutter 
oscillations (dashed line) is strongest for the integrable case 
of equal masses, whereas the first revival A 1 (solid line) grows 
for lighter impurities with a maximum at mi/m-f ~ 2/3. The 
shift of the maximum can be attributed to the formation of 
the correlation hole at short times. 



with the quantum flutter frequency. This is analogous to 
the plasmon-magnon gap discussed above, however we 
must now find the equivalent quantity for a system with 
broken translational system. 

Quantum flutter is intuitively due to beating between 
two types of states, magnon-like and plasmon-like, the 
frequency of which is given by the plasmon-magnon gap. 
In the TG regime this energy difference can be expressed 
through absolute ground state energies; quantities which 
are also meaningful for a system with broken transla- 
tional symmetry. Specifically, the magnon-like state cor- 
responds to putting the impurity in the lowest unoccu- 
pied orbital and fully symmetrizing the spin sector, thus 
eliminating all interactions. The energy of this state is 
simply Epf-f + |tz-|- since it is equivalent to the ground state 
with N+l up spin particles and a global spin rotation ap- 
plied. When the system is translationally invariant this 
state corresponds to the ground state at momentum kp. 
This can be seen by realizing that the minimum kinetic 
energy of the system, given that the total momentum is 
kp, is EN^ + kp/(2mi). The Fermi sea of N+ 1 particles, 
i.e., the state with fully symmetrized spin sector, has this 
kinetic energy, provided the masses are equal as we con- 
sider here, and no interaction energy, therefore it is the 
ground state at kp, which is the magnon-like state we 
were looking for and is a Slater determinant of the low- 
est iV + 1 orbitals. If we replace plane waves with orbitals 




FIG. 8. Harmonic confinement potential: (Color on- 
line) Quantum flutter in the presence of harmonic confine- 
ment with frequency ujt- The background gas is strongly 
interacting 7t, s = 00 and the initial impurity momentum 
Q — 1.5/cf- Other system parameters are chosen as in Fig. 4. 
The presence of the harmonic trap does not affect the quali- 
tative features, however the flutter frequency grows quadrati- 
cally with the trapping frequency, see red square in the inset. 
The gray line in the inset denotes the plasmon-magnon gap 
extracted using DMRG, see text. Thus quantum flutter fre- 
quency once again corresponds to the plasmon-magnon gap. 



that are eigenstates of the trapping potential then analo- 
gous states can be created: for example, the "Fermi sea" 
of N particles becomes the occupation of the lowest N 
eigenstates, which are Hermite functions for a harmonic 
trap, the background chemical potential becomes the en- 
ergy of the lowest unoccupied state, and the energy of a 
deep hole is the energy of the lowest occupied state. 

Meanwhile, the plasmon-like state starts with the 
background of N up spins, takes a particle from the 
lowest occupied orbital to the lowest unoccupied orbital 
making a deep hole and adds a dressed impurity in its 
ground state. The energy of the plasmon-like state con- 
sists thus of the following three contributions: (i) the en- 
ergy of the Fermi sea .Ejvt- (ii) the energy of the lowest 
unoccupied orbital fj^ minus the energy of the lowest oc- 
cupied orbital £"™ m , and (iii) the dressed impurity chem- 
ical potential /Lt7. Now, for equal masses, E™ ln is also the 

impurity chemical potential A*| 7_ °'' when the impurity 
does not interact with the background, as it will then oc- 
cupy the lowest orbital. Using this fact we obtain for the 
energy of the plasmon-like state: Epj^+fi-f— Thus 
the energy difference between plasmon- and magnon-like 
states becomes huif — fij — (J,® and can be calculated with 
DMRG, which is how we obtain the gray line in Fig. 8. 

The dependence of the flutter frequency on the trap 
frequency is well described by that energy gap. For a 
background gas consisting of Rb atoms of linear density 
p = 10 4 cm — 1 at the trap center, the simulated trapping 
frequencies correspond to top ~ 27rx{6, 18, 29} Hz (main 
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figure). 



IV. DYNAMICS IN THE PRESENCE OF A 
CONSTANT DRIVING FORCE 

Turning to the case of constant force acting on the 
impurity, we start the impurity in a plane wave with 
momentum zero and track its momentum as a function 
of time. We find that initially the impurity accelerates 
as though it was free falling and subsequently after a 
few (roughly two) Fermi times reaches a final velocity Vd 
(called drift velocity in [27]) around which it performs 
Bloch oscillations, see Fig. 9(a). These Bloch oscilla- 
tions occur due to the 2kp periodicity of the ground state 
dispersion of interacting one-dimensional systems. Thus 
they have to be distinguished from Bloch oscillations of 
free particles on a lattice, where the presence of the lat- 
tice itself renders the dispersion 2n periodic. Note that 
for infinitesimal driving, the impurity would adiabatically 
follow the ground state branch of the dispersion, and thus 
not move on average. In contrast for finite driving force 
the acceleration of the impurity is accompanied with the 
continuous emission of phonons in the background gas, 
which in turn leads to the finite drift velocity 

To quantify the acceleration G we introduce its dimen- 
sionless pendant as e a = rn 2 G/p 3 , which is related to 
the parameter F in the spin- 1/2 Bose- Hubbard model 
Eq. (1) through F — 2a 3 p 3 e a Jj-. Model parameters in 
the simulations, which are done on a dilute lattice, are 
chosen such that the period of Bloch oscillations due to 
the 27r-periodicity of the underlying single-particle lattice 
dispersion vastly exceeds the simulation time. 

We compare the frequency of Bloch oscillations ob- 
tained from numerical simulations with the semiclassical 
result Eq. (64) from Ref. [27] in Fig. 9(b). The Bloch 
oscillation frequency grows with increasing force e a as 
one would expect, since the force drives the momentum 
increase of the impurity. The results we obtain are sig- 
nificantly renormalized from the semiclassical prediction, 
suggesting higher order interaction effects are playing an 
important role. We can rule out quantum flutter as be- 
ing the culprit for the oscillations, since the energy scales 
are completely different. Note that the Bloch oscillations 
around the drift velocity are strongly damped and do 
even exist if the impurity exceeds the sound velocity of 
the background gas, see Fig. 9 (a), which is not captured 
in the semiclassical approach of [27]. Reducing the back- 
ground gas interaction from the TG regime 7^ = 00 to 
7b 9 = 5 does not significantly change the oscillation fre- 
quency, Fig. 9(b), whereas the drift velocity increases, 
Fig. 9 (a). 

For diffusive impurity dynamics the drift velocity is lin- 
ear proportional to the force Vd = fJ-F, where p is the mo- 
bility of the impurity, which is shown in Fig. 9 (c) along 
with predictions from [27]. In contrast to the semiclassi- 
cal predictions of Ref. [27], the mobility extracted from 
the numerical simulations is not independent of the force 




FIG. 9. Bloch Oscillations: (Color online) The impurity 
is accelerated by a constant force e a . After an initial free 
acceleration the impurity momentum saturates at the drift 
velocity around which it oscillates, (a) shows (Pi(t)) for 
7 = 12, background gas interactions 7^3 = 00 and 7(, 3 = 5, 
and t a = {8, 16, 24}, respectively, see legend, (b) compares 
the numerically extracted oscillation frequency uut, (c) the 
mobility p, and (d) the ratio of the drift velocity Vd and 
the oscillation frequency uii,, with the results obtained from a 
semiclassical analysis, which is valid up to e a ~ 20 [27]. 



e a . In contrast, the ratio of the drift velocity Vd and the 
Bloch oscillation frequency Wj,, which in the semiclassical 
analysis is proportional to the mobility, is approximately 
constant over a very large range of driving e a also when 
full quantum mechanical fluctuations are taken into ac- 
count (d). 

In Fig. 10(a), we compare the expected position of 
the impurity for both freely accelerated atoms 7 = 
which follow a parabolic curve and impurities interacting 
with the background gas with strength 7 = 12. Initially, 
the trajectories of the two cases cannot be distinguished. 
However, with increasing time the 7 = 12 trajectory 
departs from the initial parabolic behavior and devel- 
ops a linear trend, reflecting that the drift velocity Vd 
is reached. A renormalization of the impurity trajectory 
with interactions has also been found in experiments [1], 
where an imbalanced gas with interaction 7 = 7 has been 
created, the impurity accelerated with e a = 67, and its 
trajectory has been measured up to ~ 2tp. In Fig. 10 (b) 
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FIG. 10. Impurity trajectories: (Color online) Trajectories 
of the impurity which is driven by a constant force, (a) com- 
pares the impurity trajectories for finite interactions 7 = 12 
to free fall 7 = for different values of the acceleration e a = 4 
and 40. (b) compares without any adjustable parameters 
with the data from experiments [1], where we included the 
points that correspond to the impurity still being immersed 
in the background gas, as discussed in Fig. 3 of [1]. In the nu- 
merical simulations we fully treat both the interactions and 
the harmonic confinement. 



we compare theory and experiment without adjustable 
parameters, and find that they are quantitatively com- 
patible. Our calculations fully include both the harmonic 
trap potential, the interactions 7 and 7& ff = 00, as well 
as driving force e a . Note that the trajectory of a free 
particle in a harmonic trap follows a cosine whose pe- 
riod in the case of the experiment [1] is however much 
larger than the observed time scales, see Fig. 10 (b), yel- 
low (light shaded) curve. 

The constant force acting on the impurity gradually 
drives it along the ground state dispersion while continu- 
ously emitting phonons at the speed of sound u s [27] — a 
process which heats the background gas. The power 

?>t = >t(*)> 

dissipated into the background gas is shown in Fig. 11. 
At short times V-\ grows quickly but when the impurity 



e-a = 0,Jbg = 00 

ea=16,7b s =oo 
e Q =24, 7 &s =oo 




FIG. 11. Dissipated power: (Color online) The power Vf 
dissipated in the background gas due to the continuous emis- 
sion of phonons increases strongly at short times. When the 
impurity reaches its drift velocity Vf reaches a plateau around 
which it oscillates, providing another signature of Bloch oscil- 
lations. Times at which the impurity momentum has a local 
maximum are indicated by squares. The maxima in Vf thus 
appear at slightly later times than the maxima in the impurity 
momentum, indicating a retardation effect of the bath. 



reaches its drift velocity reaches a plateau (in clas- 
sical systems the constant average heating rate would 
correspond to Joule heating) around which it oscillates. 
These oscillations are a further manifestation of Bloch 
oscillations. The squares highlight the times when the 
maxima in the impurity momentum (P^(i)) appear. V-\ 
shows maxima at the same frequency as those in (P|_(t)), 
however they are lag behind in time. 

In addition to the heating of the background gas, the 
picture of the continuous phonon emission can be directly 
confirmed by the dynamical density-density correlation 
function G^i(x,t) [Eq. (12)], see Fig. 12, which looks at 
the density of the background liquid with respect to the 
impurity position. We find that the significant features 
occur within a spacetime cone whose apex is at the impu- 
rity position at t = 0, and which is bounded by the speed 
of the sound u s of the background gas (indicated by the 
dashed lines in Fig. 12, center column). Within the cone, 
in the forward direction, i.e., the direction in which the 
impurity is accelerated, the density is increased, while in 
the backwards direction it is decreased. Thus the impu- 
rity is constantly emitting particles and holes symmet- 
rically at the speed of sound u s around its position. In 
the TG limit 7h g = 00 the speed of sound u s is given by 
the Fermi velocity vf while for 7t g = 5 it is renormalized 
to u s ~ 0.8vf [22]. The different values of the speed of 
sound u s are directly apparent in the different angles of 
the cones. 

The dynamics of G^(x,t) reveals further details of 
the constantly accelerated impurities: Due to the re- 
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(a.l) (a.2) (a.3) 




FIG. 12. Background density for an impurity in gravity. (Color online) Time evolution of the correlation function 
G-fi(x,t) = (pi(0,t)p^(x,t)) , which represents the density of the background gas around the impurity, for 7 = 12, a finite 
external force e a = 8, and background interaction 755 = 00, upper row, and 7b s = 5, lower row. The other system parameters 
are chosen as in Fig. 4. The left column shows a perspective view on Gf±(x,t), the center column a top view, and the right 
column cuts through Gfi(x,t) at constant times. The main effect of the impurity on the background density is concentrated 
around a cone whose demarcations are the dashed lines with slopes given by the speed of sound u s ■ In the forward direction (the 
same direction as the acceleration), the density is increased, while behind the impurity the density is depleted. If we assume 
that as the impurity creates a correlation hole it emits sound wave packets in the forward and backward directions at t=0 at the 
speed of sound u s , see shortest times in the right column, then since it is moving with an expected velocity v(t) = (Pj.(t))/m_|_, 
the position of the wave packets at time t will be f Q dt'(±u s — v(t')), where the plus (minus) sign corresponds to the forward 
(backward) direction. These lines are indicated in black. In green we highlight the maximum of Gf±(x,t) for x > 0, and the 
same for x < 0. In backward direction the green points and black lines almost coincide. In the forward direction the emitted 
wave packet vanishes in the subsequently emitted flood of phonons. These signatures should be detectable experimentally when 
measuring the background density with respect to the impurity position. 



pulsive interactions, at short times a correlation hole 
is formed around the impurity and excess particles of 
the background gas are emitted almost symmetrically 
both in forward and in backward direction at the speed 
of sound u s , see t = OMp curve in the right column 
of Fig. 12. The impurity is accelerating, thus its ve- 
locity v(t) = (P±(t))/m± is changing over time. In 
its frame, therefore, the wave packets will be centered 
around J^dt'(±u s — v(t')), the plus (minus) sign cor- 
responding to the wave packet emitted in the forward 
(backward) direction. In Fig. 12 the lines following this 
argument are highlighted in black. In the backward direc- 
tion they coincide with the maximum of the emitted wave 
packets highlighted in green, where we present the case 
of a fermionic background (top row) and a background 
with finite interactions (bottom row). In the forward di- 
rection the initially emitted wave packet is overshadowed 
by the subsequent emission of phonons. These features 



could be measured in a cold atom realization by imaging 
the impurity and the background cloud, averaging over 
several realizations and measuring position with respect 
to the impurity position. 

As discussed in the previous paragraphs, in the long 
time limit t ^ 2tF the impurity reaches the drift velocity 
around which it weakly oscillates. Therefore the impurity 
frame is approximately constant provided t > 2tp, essen- 
tially allowing in that regime for the symmetric shape of 
the cone, which reflects the emitted particle-hole pairs, 
and for the approximately linear propagation of the wave 
packets created while the correlation hole is formed. 



V. CONCLUSIONS 

In summary, we have studied two types of oscillations 
arising from the nonequilibrium dynamics of an impu- 
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rity coupled to a one-dimensional quantum liquid. In 
the first case, the impurity is given an initial momentum 
from which it only sheds a part to the background liquid. 
Further the impurity momentum undergoes coherent os- 
cillations in time provided that the initial momentum is 
of the order of or larger than the Fermi momentum. In 
the second case, the impurity is initially at rest but sub- 
ject to a constant force, which accelerates the impurity 
up to a finite velocity around which it performs Bloch os- 
cillations arising from the 2kp periodicity of the ground 
state branch. The results obtained in both cases are in- 
dependent of whether the system is intcgrablc or not and 
are robust to changes in the parameters of the system, in- 
cluding a trapping potential, which lends support to the 
suggestion of seeing signatures of the physical phenom- 
ena we describe in a cold atomic setup. Our results are 
in quantitative agreement and therefore possibly an ex- 
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